Damped harmonic oscillator: pure states of the bath 
and exact master equations 
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Time evolution of a harmonic oscillator linearly coupled to a heat bath is compared for three 
classes of initial states for the bath modes - grand canonical ensemble, number states and coherent 
states. It is shown that for a wide class of number states the behavior of the oscillator is similar 
to the case of the equilibrium bath. If the bath modes are initially in coherent states, then the 
variances of the oscillator coordinate and momentum, as well as its entanglement to the bath, 
asymptotically approach the same values as for the oscillator at zero temperature and the average 
coordinate and momentum show a Brownian-like behavior. We derive an exact master equation for 
the characteristic function of the oscillator valid for arbitrary factorized initial conditions. In the 
case of the equilibrium bath this equation reduces to an equation of the Hu-Paz-Zhang type, while 
for the coherent states bath it leads to an exact stochastic master equation with a multiplicative 
noise. 

PACS numbers: 05.30.-d, 05.40.-a, 02.50.-r 
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I. INTRODUCTION 

The model of an oscillator linearly coupled to the bath 
of harmonic oscillator has played an important part in 
statistical mechanics 0, 0, Q , quantum optics 0, HJ and 
quantum measurement theory |E S ■ In most studies 
of this model the initial state for the whole system is 
taken as a mixed density matrix. In particular, one often 
uses factorized initial state where the bath modes are 
in thermal equilibrium and the oscillator is in a pure 
state. Less often the so-called thermal initial conditions 
are used 0,E3,O)0- Pure states of the bath are rarely 
considered [L| , except for the vacuum state of the bath. 

Our goal in this paper is to compare the behavior of 
the oscillator for different pure initial states for the bath 
modes to the case of the bath in equilibrium. The quanti- 
ties we will be considering are the averages and variances 
of the oscillator coordinate and momentum, as well as 
Trp 2 as a measure of oscillator entanglement to the bath. 
We also would like to show how different initial states 
for the bath modes lead to different exact master equa- 
tions for the oscillator density matrix. In deriving such 
equations we will use an exact formal solution for the 
characteristic function of the oscillator, rather then the 
path i nteg ral techniques for the reduced density matrix 
|7l fToL ITTi Il4 | . This approach makes it possible for this 
model to obtain master equations for arbitrary factorized 
initial conditions. 

This paper is organized as follows. In Sec. [H]we con- 
sider the model and its exact solution. In Sections |IIII II VI 
andlvlthe oscillator behavior is considered, respectively, 
for the bath modes in equilibrium, number states and 
coherent states. Exact master equations are discussed in 
rvll Concluding remarks are given in Sec. IVIII 



II. THE MODEL OF A LINEARLY COUPLED 
OSCILLATOR 



The system Hamiltonian is given by 
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Here and a are the creation and annihilation operators 
of the harmonic oscillator, and b\ and b k are the creation 
and annihilation operators for the bath modes. The co- 
efficients are assumed to be such that the Hamiltonian is 
a positive definite quadratic form. The coordinate and 
momentum operators for the oscillator are related to 
and a through 
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where m is the oscillator mass. By a suitable choice of 
coefficients Hamiltonian Q reduces to the Hamiltonian 
with coordinate-coordinate coupling or the rotating wave 
approximation (RWA) Hamiltonian. In particular, in the 
latter case the last two terms in JQl are dropped. We will 
assume that if the number of the bath modes increases to 
infinity the frequency io k becomes a continuous function 
of k. We will refer to the limit of the infinite number 
of modes with average energy of each mode being held 
constant as the thermodynamic limit. 

Various forms of Hamiltonian Q corresponding to dif- 
ferent choices of frequencies and coupling constants as 
well as its exact diagonalization has been extensively 
studied in the literature. General but formal discussion 
of the diagonalization of a Hermitian quadratic bosonic 
form [of which is a special case] can be found in . 
Systems of oscillators with coordinate only coupling were 
considered in 0] . Detailed investigation of Hamiltonian 
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||TJ in the case of coordinate coupling can be found in ref- 
erences 0, H, El • Relation between several forms of Q 
and the translationally invariant Hamiltonian with coor- 
dinate coupling was discussed in 

El- 

Relation between 

the coordinate coupling and RWA is studied in |l2lll7| . 

The equations of motion for the annihilation and cre- 
ation operators are 



a = -iva - i ^ Ukbk - i ^ v kb\, 

k k 

b\ = —iuJkbk — iu* k a — ivka, 
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b\ = icOkb k + iuko) + iv* k a. 
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This system of equations can be solved subject to the set 
of initial conditions a(0) = a, bk(0) = bk, a^(O) = aJ, and 
&1(0) = bt. Since system (0 is linear, its solutions will 
depend linearly on the initial conditions. In particular, 
a(t) is given by 

a(t) = A(t)a + Bk(t)b k + C{t)cJ + £ D k {t)b{. (4) 

k k 

Similarly, for a* it) we have 
ot(t) =A*(t)ai+Y,Bmbl+C*(t)a+J2Dmbk- (5) 

k k 

Coefficients A(t), Bdt), C(t), and Dk{t) satisfy the fol- 
lowing relation (see 15] for details): 
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\C(t)\ 2 +J2\B k (t)\' 2 -^2\D k (t)\ 2 = 1. (6) 



They can, in principle, be calculated for a each partic- 
ular form of frequencies and coupling parameters. For 



the purposes of this paper we will not need explicit ex- 
pressions for these coefficients. The only assumptions we 
will use is that in the thermodynamic limit coefficients 
A(t) and C(t) vanish for t — > oo, and B k {t) and Dk{t) 
remain bounded in the same limit. Physically these re- 
quirements correspond to the fact the the initial state of 
the oscillator is forgotten for long times while any ob- 
servables associated with it (e.g., average energy) remain 
finite. The detailed calculations of coefficients that show 
such behavior as well as conditions on the coupling con- 
stants and frequencies in the thermodynamic limit can 
be found in the original references 0, S 0> DU LLd] • 

The reduced dynamics of the oscillator is conveniently 
described in terms of a symmetrically ordered character- 
istic function defined as |j| 



(7) 



Here xiVi t) can be treated as a function of 77, rf and t. 
We will always suppress the second variable to simplify 
the notation. Using characteristic function J7J) we can 
calculate expectation values of the symmetrized prod- 
ucts of operators a) and a. It can also be converted to 
any of the quasi-probability distribution functions or the 
reduced density matrix We will also use Trp 2 as a 
measure of purity of the oscillator state p. In terms of 
the characteristic function Trp 2 is given by 



Trp 2 



d 2 v\x(v)\ 



(8) 



To see how xiVi *) evolves in time we use the Heisen- 
berg picture and insert expressions an( l © m t° 
to obtain 
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(9) 



The time dependence of the coefficients is suppressed here 
and whenever possible to avoid heavy notation. 

In this paper we consider only factorized initial condi- 
tions, i.e., we assume that the initial density matrix of 
the whole system ptotai factorizes into the oscillator and 
bath density matrices as ptotai = p <£> Pbath- In this case 
the characteristic function takes the form 



x(v,t) = x((vA*-v*c),o)F( v ,t), (10) 



with 

FM = T r ^JJe^J-^- *)^-^- 8 *-^!) 6 ^ . 

(11) 

Equation i|l(J[l expresses the oscillator characteristic func- 
tion at time t in terms of the initial characteristic func- 
tion. In the limit of long times, when A and C van- 
ish, characteristic function ^(77, t) is determined by the 
asymptotic form of F(t],t). 

We will now consider the time evolution and asymp- 
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totic values of the average oscillator coordinate and mo- _ B k D k ink + — ) (16) 



mentum, the coordinate and momentum variances and \ 2 

Trp 2 for different initial states of the bath 



III. EQUILIBRIUM STATE OF THE BATH 

This section consists primarily of an overview of well- 
known results. The state of the bath is the grand canon- 
ical ensemble given by 



Pe* = Hi 1 

k 



e -[3huj k \ e -phu} k b\b k 



Using the identity 



Tr [pe q e r ' 



ribl- V *b k \ _ -|r,| 2 (n fc + i) 



with 



e 0hu) k _ i > 
we find function F(j], t) to be 

F(rj,t) = e-^l'+^+i"?* 5 
Here coefficients a and 7 are given by 



a 



J2(\B k \ 2 +\D k \ 2 ){n k ' 



(12) 



(14) 



(15) 



We will often use coefficients a and 7 in the limits t — > 00 
or T — > 0. In such cases we will use the notation aioo, 700 
and a , 7 , respectively. If both limits are taken we will 
use a ^ and 7^. 

The time evolution of the average coordinate and mo- 
mentum is easily calculated either through the charac- 
teristic function, or by directly using (QJ and J5J. 



(13) (*(«)) = J^-({A* + C)(a!) + {A + C*)(a)), 

V 2mv 

(p(*)> = i\/^((A* + C)<ot) -(A + C*){a))(17) 



There is no dependence on the state of the bath. For 
arbitrary temperature (x(t)} and (p(t)} depend only on 
their initial average values. In the limit of infinitely long 
times (x(t)) and (p(t)} vanish. 

For the variances of the oscillator coordinate and mo- 
mentum we obtain 



(x 2 (t))-(x(t)) 2 = 



(A* + C) 2 ((at a t) - ( a t)( a t)) + (A + C*f((aa) - (a) (a)) 



(p 2 (t))-(p(t)) 2 = 



+2| (A + C*)| 2 «at a ) - (at)( a ) + I) + 2 (a + 7 + 7*)' 
^f[-(A* C) 2 «aV> _ (a t )(a t)) _ (A-C*) 2 ((aa) - (a)(a» 
+2| [A C*)| 2 «at a ) - (at) (a) + ~) + 2(a - 7 - 7*)^ 



(18) 



r 



These quantities depend both on the initial state of the 
oscillator and the temperature of the bath. For t — > 00 
we have 



[X ) - [x) 
(P 2 ) (P) 2 



(ttoO +7CX3 + 7»)l 

mv 

hmv(aoo - 7oo - 7^). (19) 



Coefficients a M and 700 are proportional to n k . There- 
fore, at high temperatures variances l|19H grow as kT. 

The measure of the oscillator purity is given by the 



integral: 

Trp 2 = i y'd 2 7 ? |x((77^*-r7*C , ),0)| 2 e- 2 ( Q l"l 2 - 7 * 1 ' 2 - 7 ''* 2 ). 

(20) 

This integral can be calculated for specific initial states of 
the oscillator. For infinitely long times, when the initial 
state is forgotten, we obtain 

1 



Trp 1 



2 A /a2 -4| 7 oo| s 



(21) 



For high temperatures Trp 2 is proportional to 1/kT . The 
state of the oscillator becomes less pure as the tempera- 
ture grows. Let us note that for intermediate times Trp 2 
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can take lower values than its value at t — * oo. For zero 
temperature of the bath the asymptotic value of Trp 2 
will, in general, be less then one since the oscillator re- 
mains dressed at T = 0. 

For the special case of the RWA Hamiltonian, coeffi- 
cients C(t) and Dk(i) in and JSJ) are equal to zero. As 
a result, = | and 7 = and for long times Trp 2 = 1. 
In this case the reduced vacuum of the oscillator is a pure 
state that is identical to the ground state of the uncou- 
pled oscillator. 



IV. THE NUMBER STATES FOR THE BATH 
MODES 

We now consider the initial state of the bath with each 
mode in a number state. 

|{rife}> = \n kl )® |nfc 3 )<g>--- (22) 

with {n^} denoting a set of occupation numbers n k for 
all modes. 



Function F(rj. t) can be calculated using the identity 


(n k \e^-^\n k )=e-^L nk (\ V \ 2 ), (23) 

where L n (x) is a Laguerre polynomial. We obtain 

F{ V ,t) = e -«°M 2 +7°V+7 l V 2 

xnMl(»?B**-rfD fc )| 2 ). (24) 

k 

The behavior of average coordinate and momentum of 
the oscillator is exactly the same as for the equilibrium 
state of the bath and given by (|17|l . Thus, (x(t)) and 
(p(t)) do not depend on the particular set of occupation 
numbers. Both quantities vanish for t — > 00. For the 
variances of the coordinate and momentum we have 



(x 2 (t))-( x (t)y 



Ivnv 



(A* + C) 2 ((aV) - (a t )(a t )) + (A + C*f({aa) - (a)(a)) 



+2| (A + C*)\ 2 {(a^a) - (<J)(a) + \) + 2(5 + 7 + 7*) 
hmv 



(p 2 (t))-{p(t)y = — -(A* - C) ((a t a t ) - <a t )(a t )) - (A-C*) ((aa) - (a) (a)) 



+2\(A - C*)| 2 ((at a ) - < a t)( a ) + I) + 2(5 -7 - 7*) 



r 



(25) 



where 



= J2{\ B k\ +\ D k 

k 



n k + - 



(26) 



The variances differ from the equilibrium ensemble case 
by the replacement of njt with rik ■ 

To get a better picture of how these values relate to the 
equilibrium case we have to make some assumption about 
the occupation numbers. Let us consider an ensemble 
of the occupation number states corresponding to the 
number state decomposition of the equilibrium density 
matrix 



(27) 



{«fc} 



Here the probability for a particular set of occupation 



numbers P n ({nk}) is given by 



Pn({n k }) =I1( 1 - 



e -0hui k \ e -/3huj k n k 



(28) 



We now assume that the number states arc taken from 
ensemble (12 7|) . Any quantities calculated for each indi- 
vidual number state (e.g., averages, variances, Trp 2 ) can 
then be averaged over P n {{nk}) to obtain their average 
values in ensemble (|27|l . These latter averages will give 
typical values for the pure state quantities in ensemble 

G3- 

Averaging variances H25|) over ensemble (|27|l will give 
the same variances as for the equilibrium case. In par- 
ticular, in the limit of long times we have 



(x 2 ) - (x) = (««, +7oo +7^)> 

(p 2 ) - (p) 2 = hmviaoo - 700 - 7^). 



(29) 



Here we use overlining to denote averaging over the en- 
semble of pure states. 
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The fact that the average variances are the same as for 
the equilibrium ensemble is not surprising since expres- 
sions in 125(1 are linear in n k 's. More importantly, in the 
thermodynamic limit almost all states in ensemble 1(27(1 
will have the same variances as in equilibrium. Let us 
show this for the coordinate variance. We can treat the 
coordinate variance as a function of random variables n k 
described by the distribution P n ({n-k}) ■ Calculating the 
variance of this function for the distribution P n ({n k }) 
we obtain. 

((x 2 ) - (x) 2 ) 2 - ((x 2 ) - (x) 2 ) ■ ((x 2 ) - (x) 2 ) 

= ^EK^+^)I 4 (^-^ 2 )- (30) 

k 



If number states are taken from the ensemble ((27(1 we can 
calculate average Trp 2 for number states in this ensem- 
ble. Using the identity [2(j 

E^^Y^exp^^l), ( N<1 ), 

n 

(32) 

we obtain 

Trp 2 = ~ y'd 2 7 ? |x(M*-77*C),0)| 2 e- 2 ( tt l' ) l 2 -^*'' 2 -^* 2 ). 

(33) 

This is exactly the same as Trp 2 in ((20|) . We can see 
that, at least on the average, Trp 2 for number states 
from ensemble ((2711 is the same at all times as for the 
case of equilibrium bath. As a consequence, in the limit 
of long times the state of the oscillator becomes less pure 
for higher occupation numbers for the bath modes. 

V. COHERENT STATES FOR THE BATH 
MODES 

We now consider the case where all bath modes are 
initially in coherent states. 

\{flk}) = \Pki)®\Pk a )®--- (34) 

with {f3k\ denoting a set of complex numbers /3& speci- 
fying the coherent states. One can interpret such a state 
as the most classical state of the bath. Function F(r), t) 
is calculated to be 

F(rj,t) = e <5>-<V-« l'!l 2 +7°V+7V 2 . (35) 

Here 

5 = Y J {BA + D k p* h ). (36) 

k 



We note that coefficients B k and D k must depend on 
the number of bath modes iV as 1/ y/N in order for the 
quantities like 1(25(1 to remain finite in the thermodynamic 
limit. Therefore the sum over k in l(3U(l is proportional to 
1/./V and vanishes for N — > oo. A similar argument can 
be applied to ((p 2 (t)} - (p(t)) 2 ). Thus, as in equilibrium 
case, we expect the variances to grow as kT for high 
occupation numbers. 



Let us now consider the behavior of Trp 2 . Using defi- 
nition ||SJ and characteristic function l(24|) we obtain 



(31) 

I 

For the average coordinate and momentum we obtain 

<*(*)) = \[^((A*+C)(al) + (A + C*)(a) 
+6 + 5*), 

(p(t)) = is J^({A*-C)(<J)-{A-C*)(a) 

+5* -5). (37) 

To get a better understanding of how (x(t)) and (p(t)) 
behave we have to make some assumption about param- 
eters j3 k . 

Analogously to the procedure used for the number 
state bath let us assume that the coherent states are 
taken from the ensemble corresponding to the coherent 
states decomposition of the equilibrium density matrix 

p eq = J d 2 {(3 k }\{(3 k })P c ({(3 k })({(3 k }\. (38) 
The probability distribution P c ({(3k}) is given by 

P c ({/3 fc })=n^ ex pf Z ^V ( 39 ) 
A A irn k \ n k J 

Decomposition (13811 is just the P-representation of the 
equilibrium density matrix 0]. 

If the coherent states are taken from ensemble ((38II . 
then, for each set of the coherent states, 8{t) is a realiza- 
tion of a complex colored normal noise with zero mean 
and correlation functions given by 

SWW) = E (Mt)Bt(t') + D k (t)Dt(t'))n k , 

k 

WO = E (Bk(t)Dk(t') + B k (t')D k (t)) n k . (40) 

k 



I 

Trp 2 = I / ^^((^-^C)^)^" 2 ^!''! 2 ^ *, 2 ^ 2 )^^ (K^.^f) . 
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Clearly, (x(t)) and (p(t)) are also realizations of the nor- 
mal noise. The mean in this case is, in general, non-zero 
but will go to zero for long times. 

Let us evaluate the typical asymptotic values taken by 
(x(t)} and (p{t)). This can be done by calculating average 
(x) and (p) for ensemble (13811 . 



(pY 



"00+700 + 7* 



Timv ( a 



loo - 7o 



7° ' 

I OO 



^*oo 1 / OO 1 I c 



(41) 



(42) 



For high temperatures these quantities are proportional 
to kT . Therefore typical asymptotic values for (x) and 
(p) are proportional to \/kT. 

The coordinate and momentum variances are the same 
as for the case of equilibrium bath at zero temperature, 
i.e., they are given by l|18fl with a and 7 replaced by a 
and 7 , respectively. Therefore, these variances do not 
depend on a particular set of parameters (3k and their 
asymptotic values are the same as for the reduced vacuum 
state of the oscillator. 

The value of Trp 2 is given by 



Trp A = 



2 = - fd 2 V \ X ((vA*-ifC),0)\ 2 e^ aoM2 -^ 2 -iov' 2 )^ 



(43) 



and is the same as in the case of zero temperature bath 
with the asymptotic value given by 



Trp 2 



2V(ag o r-4hg 



(44) 



More generally, the asymptotic characteristic function 
is determined by the asymptotic form of function F(rj, t) 
(|35[) and corresponds to the density matrix that is a dis- 
placed reduced vacuum state p vac , 

p = D(S) Pvac DHS), (45) 

where the displacement operator D(8) is given by 

D(S) = e Sai - s ' a . (46) 

We can conclude that for coherent states bath in the 
limit of long times the oscillator becomes localized in 
the following sense: the coordinate and momentum vari- 
ances, as well as entanglement to the bath, are same as 
for the reduced vacuum (and do not depend on tempera- 
ture), while the average coordinate and momentum ran- 
domly fluctuate with the amplitude of the fluctuations 
proportional to \/kT. The rate of such localization is of 
the order of the relaxation rate as can be concluded from 
(QUI for T = 0. 

For the special case of the RWA Hamiltonian the 
asymptotic state of the oscillator is a pure coherent state. 
The oscillator can start in a mixed or pure state, it then 
goes through a period of entanglement and asymptoti- 
cally its state becomes a pure coherent state with ran- 
domly fluctuating average coordinate and momentum. If 
the oscillator is initially in a coherent state it will al- 
ways remain in a coherent state, never entangling with 
the bath |2l| . In Appendix we give explicit expression 
for Trp 2 (t) for two types of initial oscillator states for 
the RWA case. 



VI. EXACT MASTER EQUATIONS 

There is a similarity of a mathematical nature between 
the coherent states bath and bath in equilibrium. In both 
cases function F(rj,t) is Gaussian, and, as a result, the 
characteristic function satisfies simple master equations 
as will be shown shortly. 

Before doing so let us note that for any factorized ini- 
tial conditions (for any system-bath model) there always 
exists an exact equation for the reduced density matrix 
with the time evolution governed by a time-dependent 
operator that does not depend on the state of the system. 
Obtaining an explicit expression for such an operator can, 
in general, be a difficult task. Let us show that such an 
operator can be constructed for the present model. We 
will continue to use the characteristic function space be- 
cause of its mathematical convenience. If necessary the 
equations can be transformed into other representations. 
To simplify the demonstration we introduce the following 
transformed characteristic function 



(47) 



Using expression 110(1 for x(Vi t) and the fact that 
x(t7,0) = x(??,0) we can write function x(?7,i) at time 
t in terms of the x(?7, 0) as 



(48) 



We now use the fact that xi^i t) depends on 77, 77* and t 
only through (r]A* - 77* C) and (rj* A - rjC*). Differenti- 
ating solution (|48|l with respect to time we have 



dx 

dt ~ d(r]A* - r;*C) 

dx 



+ 



d{rj*A - 77C*) 



(tjA* - rfC) 
{r 1 *A-r 1 C*). (49) 
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For derivatives with respect to rj and 77* we obtain Here and are given by 

= dx v d X » 

dr? d(rjA*-Tj*C) d{r)*A-r)C*) ' 

a??* <9(rM* -77*c) 0(nM-nC*) ■ y > . . , . 

rr - — <rn - ca-ac 

Solving the last two equations for the derivatives of x \A\ 2 — \C\ 2 ' |A| 2 — |C| 2 ' 

with respect to (7/^4* — ifC) and (77* A — r/C*) and sub- 
stituting them in 14911 we obtain a closed equation for 
X- 

dx dx dx 

~di = ^ ^ + ' + + ftf • Substituting definition of X 63 into equation JHH we 

(51) obtain the following equation for x 



dx 
dt 



dlnF 

dt 



dx 
di] 



dliiF 
di] 



x) + mv*+c(t)v) 



dx 
drj* 



dlnF 

drj* 



(53) 



This is a closed equation for characteristic function 
x(i],t). The explicit form of the time-dependent oper- 
ator is determined by function F(rj,t), which, in turn, 
is determined by the initial state of the bath. We will 
not go into the analysis of this equation and limit our- 
selves to a few remarks. This equation is always local 
in the 77-space. However, transforming it into the quasi- 
probability distribution space or coordinate or momen- 
tum representation for the density matrix will, in general, 
lead to non-local equations. Only for some special forms 
of F(r), t) can we expect to get local equations in these 
representations. Note that exact master equations for 
the oscillator-bath model that appeared in the literature 
use equilibrium initial state for the bath 0, j squeezed 
equilibrium for the bath [T^j , or modified equilibrium for 
the whole system 0,^]. When converted into the char- 
acteristic function space all these equations contain time 
dependent operators that are at most quadratic in 77, 7/*, 
d/dr], and djdrf . This form of the time dependent op- 
erator is related to the fact that propagating function for 
such initial conditions is Gaussian. Inspection of equa- 
tion l|53[) shows that such a simple form of the time de- 
pendent operator in our case is possible only for Gaussian 
functions F(r),t). 

We now consider particular form of equation (|53|l for 
this special case, i.e., when function F(r),t) is given by 



F{T),t) 



a 6*r]— Sri* -a\r]\ 2 +y* r] 2 +yr]* 



(54) 



Here 5, a and 7 are time dependent parameters char- 
acterizing each particular Gaussian state of the bath. 
Function 1541) will include such states for the bath modes 
as equilibrium, squeezed states, squeezed and displaced 



equilibrium, etc. Substituting (|54|l into (|53|1 we obtain 
after some simplifications 

+ K (t)\ v \ 2 x + mWx + »(t)v* 2 x 

+**(t)r} X -<T(t)ri*X- (55) 
Coefficients n(t), n{t) and a(t) are given by 

«(*) = «(£ + 0-2(C7* + C*7)-a 

a(t) = Of -&+S. (56) 

In the case of the equilibrium bath we need to put a = 
a, 7 = 7 and <5 = in (|56|l . In this case the last two terms 
on the right-hand side of equation (|5*5|l disappear and 
the equation becomes the characteristic function version 
of the Hu-Paz-Zhang equation 0, f° r the Hamiltonian 
0, i- e., it will allow for the possibility of momentum- 
momentum and momentum-coordinate coupling between 
the bath and the oscillator. 

If the bath modes are initially in coherent states then 
we have to use a = a , 7 = 7° and 5 = 8. If the coher- 
ent states are taken from ensemble Q38fl. then, for each 
particular set of states, a(t) will be a realization of a 
complex colored normal noise. The equation becomes a 
stochastic master equation with a multiplicative noise. 
In this case (|55|) describes the oscillator localization in 
the sense discussed above with the average coordinate 
and momentum subject to random fluctuation. Equa- 
tions of such type can be of interest in the theory of 



8 



quantum measurement as an alternative to the stochas- 
tic Schrodinger equations j2^. In particular, compared 
to the latter equations, equation Ij55(l does not conserve 
purity of the state. Such behavior is more physically 
plausible for a system coupled to a bath. 



VII. CONCLUDING REMARKS 

Equilibrium density matrix is often considered as an 
irreducible concept, viz., it is assumed that this is a true 
state of the bath in each individual experiment. We be- 
lieve that such a view is an oversimplification, and it is 
more realistic to assume that true state of the bath is a 
density matrix or even a pure state which includes some 
random component. Ensemble of such states for different 
realizations of the random component gives an equilib- 
rium density matrix. Such a view is partially supported 
by classical statistical mechanics. Indeed, most physicists 
agree that the "true state" of the classical bath in thermal 
equilibrium is a point in phase space. The coarse-grained 
macroscopic observables for such a point (e.g., number of 
particles in a volume element) are essentially the same as 
calculated for one of the Gibbs ensembles (of such points) 
for large enough course-graining pflj . 

It is true that the ensemble decomposition of the equi- 
librium density matrix in terms of other density matri- 
ces (or pure states) is not unique, and that the averages 
calculated for any decomposition are identical to equi- 
librium averages. We believe, however, that it can be 
of importance, that, when quantities like variances and 
Trp 2 are calculated for individual members of the en- 
semble, we can get drastically different behavior for the 
system for different decompositions, as was shown in this 
paper. 
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APPENDIX 

In this Appendix we give explicit expressions for 
Trp 2 (t) for two types of initial pure states of the os- 
cillator for the RWA Hamiltonian if the bath modes are 
initially in coherent states. 

Let the oscillator be initially in a squeezed state ob- 
tained by acting on the oscillator ground state with the 
squeezing and displacement operators 



\ 1 p)=D(a)S(e)\0)._ 
Here the squeezing operator S(e) is given by 



S(s) 



e 2 



(A.l) 



(A.2) 



and e = re 2 " 1 ^ is a complex squeezing parameter. Using 
expression (143(1 for Trp 2 we obtain 



Trp 2 (t) = 



yJl + 4\A(t)\ 2 (l- \A(t)\ 2 ) sinh 2 r 



(A.3) 



The maximum entanglement [or minimum Trp 2 (i)] is 
reached when |A(i)| 2 = 1/2. Moreover, Trp 2 (t) is a 
symmetric function of |A(i)| 2 with respect to this point. 
When \ A(t)\ 2 = 1/2 we have 



Trp 2 = 



yl + sinh 2 i 



(A.4) 



We have larger entanglement for larger squeezing with 
Trp 2 going to zero when r goes to infinity. 

We now consider an initial state of the oscillator given 
by a superposition of two coherent states 



1 



N 



(A.5) 



where N is the normalization constant 
AT = 2 + (a\f3) + (/3\a) 
In this case we obtain for Trp 2 (t) 
2 

W 

\a-p\ A _ j\ 



(A.6) 



Tr P 2 (t) = i + -^( e -(i-w*)i 2 )i«-/3i 2 + e -\m\ 2 \<*-e\ 



(A.7) 



One can easily verify that Trp 2 {t) takes its minimum 
value again at |A(t)| 2 = 1/2, and Trp 2 {t) is again a 
symmetric function of |yl(£)| 2 with respect to that point, 
which may be a general property of this model. Let us 
consider Trp 2 {t) at \A(t)\ 2 = 1/2. We have 



Trp = 1 - — [ 1 — e 



(A.8) 



Using (|A.6|I . writing (a\(3) as Re llp and remembering that 



e 2 



we obtain 



iv = 1- (1 -* )2 



2(l + i?cos<p) 2 



(A.9) 



(A.10) 



We can see from the last expression that the minimum 
Trp 2 decreases with the decreasing overlap between the 
two coherent states. Note, however, that in this case 
Trp 2 can never become less then 1/2. 
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